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ABSTRACT 


Effects  of  surface  heat  transfer  on  boundary-layer  transition 
are  analyzed  in  a  three-part  study  using  the  Saf fman-Wilcox 
transition  model.  In  the  first  part  of  the  study,  model  pre¬ 
dictions  are  compared  with  experimental  data  for  cooled  and 
heated  aerodynamic  boundary  layers  on  smooth  flat  surfaces  and 
for  cooled  aerodynamic  boundary  layers  near  the  stagnation 
point  of  a  roughened  blunt  body.  Consistent  with  measure¬ 
ments,  tne  model  predicts,  on  the  one  hand,  that  heating 
destabilizes  a  smooth-surface  aerodynamic  boundary  layer 
and,  on  the  other  hand,  that  cooling  destabilizes  a  rough- 
surface  aerodynamic  boundary  layer.  Differences  between 
predicted  and  measured  transition-point  locations  are  with¬ 
in  experimental  error  bounds.  Then,  incipient  transition 
conditions  arc  determined  for  a  small,  heated  hydrodynamic 
body.  Again  model  predictions  agree  with  measurements  which 
indicate  that  relatively  small  amounts  of  surface  heating 
have  a  strong  stabilizing  effect  on  hydrodynamic  boundary 
layers.  In  the  final  part  of  the  study,  transition  location 
is  determined  for  a  large  hydrodynamic  body;  results  indi¬ 
cate  that  large  surface  heating  rates  are  not  substantially 
more  effective  than  smaller  rates.,  i 
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NOTATION 


SYMBOL 


DEFINITION 
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Pressure  coefficient,  (p-p^J/C ^/^pU^2 ) 

Specific  heat  at  constant  pressure 
Specific  turbulent  energy 

°arsons-Goodson— Golds chmied  shape  parameter 
Nondimensional  self-similar  velocity  profile 
0  for  two-dimensional  flow;  1  for  axisymmetric  flow 
Roughness  height 

Parsons-Goodson-Goldschmied  shape  parameter 
Pressure 

Laminar  and  turbulent  Prandtl  numbers 
Surface  heating  rate;  average  value  of  qw 
Radial  distance  from  symmetry  axis 
Parson-Goodson-Goldschmied  shape  parameters 
Empirical  constant 

Reynolds  number  based  on  arclength,  plate  length 

Reynolds  number  based  on  momentum  thickness 
Incipient  transition  Reynolds  number  based  on  V 
Arclength 

Parsons-Goodson-Goldschmied  shape  parameter 
Parsons-Goodson-Goldschmied  shape  parameter 
Temperature 

Freestream  total  temperature 
Freestream  turbulence  intensity 
Velocity  component  in  x,y  direction 
Boundary-layer-edge  velocity 
Freestream  velocity 
Body  volume 
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Distance  parallel  to,  normal  to  body  surface 
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Empirical  parameters 

Empirical  constants 
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1.  INTRODUCTION 


In  recent  years,  substantial  progress  has  been  made  toward 
maintaining  laminar  flow,  and  hence  low  drag,  over  aerodynamic/ 
hydrodynamic  bodies  at  practical  flow  speeds.  In  the  case 
of  aerodynamic  bodies,  the  use  of  carefully  designed  air¬ 
foil  shapes  (e.g.  to  provide  strong  favorable  pressure  gra¬ 
dient)  and  the  use  of  boundary-layer-control  techniques 
(e.g.  suction)  has  successfully  delayed  transition  to  chord- 
length  Reynolds  numbers,  ReSt,  of  the  order  of  10  million. 

Even  larger  values  of  Res^_  have  been  achieved  for  small  hydro- 
dynamic  bodies  through  the  use  of  surface  heating  and  pres¬ 
sure  gradient.  For  hydrodynamic  bodies,  extrapolations  based 
on  linear  stability  theory  indicate  that,  with  practicable 
amounts  of  surface  heating,  values  of  ReSt  in  excess  of  200 
million  may  be  possible  on  relatively  large  hydrodynamic 
bodies1 . 

The  reduction  in  drag  which  can  be  achieved  by  maintaining 
laminar  flow  over  any  vehicle  is  attractive  because  of  the 
reduced  power  requirements  to  move  the  vehicle.  However,  a 
penalty  is  generally  paid  in  maintaining  laminar  flow  on 
lifting  bodies  in  that  a  laminar  boundary  layer  separates 
much  more  easily  than  does  a  turbulent  boundary  layer,  re¬ 
sulting  in  a  significant  reduction  in  lift.  Hence,  main¬ 
taining  laminar  flow  is  practical  mainly  for  nonlifting 
bodies.  Submarines  and  torpedos  or,  more  generally,  hydro- 
dynamic  bodies  fall  into  the  latter  class. 

This  study  focuses  on  the  observed  pronounced  effects  of 
surface  heat  transfer  on  boundary  layer  transition.  Most 
importantly,  this  report  includes  transition  predictions 
based  on  a  relatively  new  transition  theory  for  small  and 
large  heated  hydrodynamic  bodies. 
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Section  2  presents  the  transition  equations,  including  a 
modification  needed  to  improve  transition-prediction  ac¬ 
curacy  when  surface  heat  transfer  is  present.  Included  in 
Section  3  are  transition  computations  for  various  aero¬ 
dynamic  and  hydrodynamic  boundary  layers.  The  concluding 
section  summarizes  results  and  conclusions. 


2.  FORMULATION 


2.1  EQUATIONS  OF  MOTION 

The  Saffman-Wilcox  transition  model2  ^  is  the  basic  tool 
used  in  this  study  to  analyze  effects  of  surface  heat 
transfer  on  boundary-layer  transition.  The  model's  accuracy 
has  previously  been  demonstrated  for  a  wide  variety  of  flows 
ranging  from  incompressible  boundary  layers  to  hypersonic 
blunt-body  flows.  These  applications  have  tested  the  model's 
ability  to  predict  transition  sensitivity  to  effects  of  free- 
stream  turbulence,  suction,  surface  roughness  and  pressure 
gradient.  In  all  cases,  accuracy  acceptable  for  most  en¬ 
gineering  purposes  has  been  obtained. 

For  incompressible  boundary  layers  (l.e.,  for  very  small  Mach 
numbers),  the  equations  of  motion  which  constitute  the  tran- 

sition  model  are: 

Mass  Conservation 


H<rJu>  +  W(r,Iv> 


(1) 


Momentum  Conservation 


3u  .  3u  _ 

ou37  +  pv  37  ‘ 


d£  ,  1  JL  [*rJ  (u  +  -)  —1 

'  -dS  +  J  9y  L  n  9yJ 


(2) 


Energy  Conservation 
3 (C  T)  3(C  T) 

+  pv-5r" 


3 


Turbulent  Energy 


3e  + 

pu^  + 


pv 


0  =  [e.*l0|-6*Pfl]pe  +  ^^[r3(lJ  +  a*l)3y] 


(U) 


Turbulent  Dissipation  Rate 


pu 


II1  +  pv  =  [cxlf^l  -  6pn]p£22  -  +  ] 


(5) 


m  Equations  1-5,  x  and  y  are  orthogonal  coordinates  parallel 
to  and  normal  to  a  body  surface.  The  quantity  r  is  e  ra 
coordinate  from  the  body's  symmetry  axis  while  .1=0  for  wo 
dimensional  flow  and  J=1  for  axisymmetric  flow.  The  veloc 
components  In  the  x  and  y  directions  are  denoted  by  u  an  v, 
fluid  density,  temperature,  pressure,  specific  heat  an, 
viscosity  are  denoted  by  p,  T,  p,  Cp,  and  u  respec 
quantities  e  and  0  are  specific  turbulent  energy  and  .peomc 
turbulent  dissipation  rate;  their  ratio,  e/Sl,  is  e  ur 
eddy  viscosity. 

Both  air  and  water  are  considered  in  this  study.  The  various 
thermodynamic  properties  for  air  are  related  through  the  per¬ 
fect  gas  law  and  the  Sutherland  viscosity  law.  ppen 
lists  the  pertinent  thermodynamic  properties  of  water. 

Several  empirical  parameters  appear  in  Equations  1-5.  For 
transitional  flows,  past  studies2'4  have  established  the 

following  values: 


PrT  *  0.89 


6  =  0.15 
a  =  0 . 50  » 


,  B*  =  0.09 
a*  =  0.50 


(6) 


(7) 


a*  =  0.30  [1  -  (1  -  X* )  (1  -  ReT/R0)  *'  "  ReT/Ro)-' 

a=  .1638  [1  -  (1  -X)(l  -  ReT/RQ)H(l  -  ReT/RQ)]  (8) 

where  H(x)  is  the  Heaviside  step function  and  Re^  is  turbulent 
Reynolds  number  defined  by 

Re^  =  e/fiy  (9) 

i 

Finally,  the  parameters  A#  and  Rq  are  constants  whose  values 
are 


a*  =  0.105  ,  :>0  =  0.10 


(10) 


The  quantity  A  depends  upon  the  freestream  turbulence  level, 
T* ,  and  pressure  gradient  parameter,  A,  defined  by 


T' 


and 


^e  dj_  dU. 
Pw  ve  dx 


(11) 


(12) 


The  functional  dependence  of  A  upon  A  and  T'  is 

A=  .  I05(l  +  2H(A)  j  1  -  exp  [ --40  A  ]  }  exp  [-3T'2])  ^3) 

Equations  1-3  are  the  time-averaged  conservation  equations 
with  classical  eddy-viscosity  and  eddy-heat-dif f usivity  clo¬ 
sure  approximations.  Equations  4  and  5  are  -onlinear  diffusion 
equations  for  e  and  JJ  which  provide  a  description  of  the  growth 
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of  disturbances  in  a  laminar  flow  up  to,  through,  and  beyond 

2 

transition.  Wilcox  presents  a  thorough  discussion  of  the  way 

in  which  the  model  equations  are  used  to  predict  boundary- 

layer  transition.  Boundary  conditions  suitable  for  flow  over 

both  smooth  and  rough  surfaces  are  given  by  Wilcox  and 
h 

Chambers  . 


2.2  NEUTRAL  STABILITY  CONSIDERATIONS 


Qualitative  features  of  model-predicted  transition  can  be 
conveniently  determined  by  dropping  the  convection  and  dif¬ 
fusion  terms  in  the  turbulent-energy  equation  (Equation  4 ) . 
In  doing  this  we  can  determine  the  neutral-stability  point 
which  is  defined  as  the  point  in  a  boundary  layer  where  tur¬ 
bulent  energy  generation,  ot#pe  |  3u/9y  |  ,  just  balances  turbu¬ 
lent  energy  dissipation,  8*p2fte.  Hence,  neutral  stability 
is  defined  by  the  following  condition: 


m?x 


a* 

B* 


9u/3y j  n  1 

pft  1 


(1*0 


Then,  noting  for  laminar  boundary  layers  that  ft  and  a*  are 
approximately 


ft 


20  y 
6  P^y2 


(15) 


ct#=  0 . 30  A# 


Equation  14  simplifies  to 


mf  X 


y 2  | 9u |  _  40 
v  197!  ~  ITT 


(lo) 


(17) 
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For  zero-pressure-gradient  boundary  layers  in  air,  the  velocity 
distribution  can  be  written  as 


u  =  UeF(n)  (18) 


Using  Equations  18  and  19,  Equation  17  simplifies  to  the 
following: 

(20) 

Now  for  incompressible  flow,  F(n)  is  only  weakly  dependent 
upon  surface  heat  transfer.  Hence,  Equation  20  implies  that 
the  neutral  stability  Reyi  olds  number,  Uex/vw,  varies  inversely 
as  A*2.  Experiments5  indicate  that  transition  Reynolds  number, 

p 

U  x,/v  ,  is  inversely  proportional  to  (T  /T  )  .  Assuming  that 
e  u  w  w  tr 

the  transition  Reynolds  number  is  proportional  to  the  neutral 
stability  Reynolds  number  (as  it  often  is)  implies  that  A* 
should  be  proportional  to  Tw/Tp . 

Computations  with  A*  independent  of  Tw/T£  verified  that  such  a 
dependence  is  needed.  The  unmodified  model  predicted  that 
heating  (cooling)  has  a  stabilizing  (destabilizing)  effect  on 
aerodynamic  boundary  layers,  in  contrast  to  the  measured  de¬ 
stabilizing  (.stabilizing)  effect.  These  predictions  were  un- 
surprising  as  Shamroth  and  McDonald^  find  a  similar  reversal 
in  predicted  transition  sensitivity  to  heating  and  cooling 
with  their  turbulence-model  transition  method.  Shamroth  and 
McDonald  resolve  the  problem  by  making  a  parameter  similar  to 


X*  an  increasing  function  of  T/T  .  While  this  model  re- 

W  0 

vision  is  sufficient  to  yield  accurate  transition  location 
for  aerodynamic  boundary  layers,  it  is  inadequate  for  the 

hydrodynamic  case.  That  is,  according  to  Equation  20,  tran- 

2 

sition  location,  x.  ,  is  proportional  to  y  /X#  .  For  a 

b  _  0  W 

liquid,  y  varies  as  T  “  near  room  temperature  (see  Appen- 

dix  « .  ..ence,  assuming  X*  is  proportional  to  T..  would  ir.rlv 

that  x  varies  as  T  Therefore,  increasing  surface  ten- 

t  w 

perature  would  decrease  ,  i.e.,  heating  would  destabilize 
a  hydrodynamic  boundary  layer.  Since  heating  stabilizes 
such  a  boundary  layer,  a  different  modification  is  clearly 
needed . 

The  parameter  X*  should  more  appropriately  depend  upon  y^/y^. 
On  the  one  hand,  y  /y  increases  as  T  /T  increases  for  air, 

W  G  W  0 

while,  on  the  other  hand,  y,r/yQ  decreases  as  T  /T  increases 
for  water.  Thus,  if  X*  were  proportional  to  some  power  of 
yw//ye’  the  mode-1-  could  accurately  predict  effects  of  sur¬ 
face  heat  transfer  on  both  aerodynamic  and  hydrodynamic 

o/h 

boundary  layer  transition.  Since  y^T  for  air,  the  fol¬ 
lowing  revised  form  of  Equation  .0  is  proposed: 

x*  =  0.105  (yw/yp)i4/3  (10a) 

An  additional,  less  obvious,  modification  is  needed.  Note 
that  the  parameter  X  partially  control®  the  rate  at  which 
disturbances  are  amplified  beyond  the  neutral  stability 
point.  Furthermore,  if  X  ever  becomes  sufficiently  large 
(relative  to  X*)  so  that  dissipation-rate  production, 
apft2 | 9u/9y | ,  overtakes  dissipation-rate  dissipation,  Bp2^3, 
(see  Equation  5)  before  the  neutral  stability  point  is 
reached,  transition  may  never  occur.  Thus,  if  X  remains 


8 


unaltered,  the  possibility  exists  that  a  finite  amount  of 
cooling  (heating)  will  cause  an  aerodynamic  (hydrodynamic) 
boundary  layer  to  forever  remain  laminar,  a  physically  un¬ 
realistic  prediction.  Hence,  Equation  13  must  be  repxaced 
by: 

X=.  105(yw/yp)i</3  {  1  +  2H(A)[1  -  exp(-40A)  lexp(-3T’ 2 )}  (13a) 

All  of  the  computations  presented  in  the  next  section  hav^ 
been  performed  using  the  transition  model  defined  by  Equations 
1  through  9,  10a,  11,  12,  and  13a;  as  noted  earlier,  appro¬ 
priate  surface  boundary  conditions  are  given  by  Wilcox  and 

Chambers1* . 
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3.  APPLICATIONS 


In  the  first  part  of  this  section,  the  transition  model  is 
used  to  predict  transition  sensitivity  to  surface  heat  trans¬ 
fer  for  smooth-  and  rough-surface  aerodynamic  boundary  layers; 
computed  transition-point  locations  are  compared  with  cor¬ 
responding  experimental  data.  Then,  effects  of  surface 
heating  on  a  small  hydrodynamic  body  are  computed;  qualita¬ 
tive  comparisons  are  made  with  experimental  data.  Finally, 
transition  location  is  predicted  on  a  large  hydrodynamic 
body . 

HEATED  AND  COOLED  AERODYNAMIC  BOUNDARY  LAYER^. 

One  of  the  easiest  of  all  flows  to  analyze  is  the  incompres¬ 
sible  flat-plate  boundary  layer  (FPBL) .  Furthermore,  analyz¬ 
ing  this  flow  provides  a  good  test  of  the  transition  model 
as  detailed  measurements  have  been  made  to  determine  transi¬ 
tion  sensitivity  to  surface  heat  transfer.  Using  an  incompres¬ 
sible  version  of  DCW  Industries'  EDDYBL  computer  code  ,  tran¬ 
sition  computations  were  performed  for  an  incompressible  PPBL 

with 

0.5  <  T/T  <  3.0  (21) 

As  in  all  computations  in  this  study,  the  value  of  ft  at  the 
boundary-layer  edge,  ,  was  given  by 

=  .Ol85U*/ye  (22) 

a  value  generally  used  in  EDDYBL  transition  calculations.  Re¬ 
sults  of  the  computations  are  shown  in  Figure  1;  experimental 
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data  of  Zysina  -  Molozhen  and  Kuznetsova  are  included  in 
the  figure  for  comparison. 

Because  the  freestream  turbulence  level  in  the  experiments  was 
not  available,  computations  were  done  for  a  relatively  high 
level  (T'=l. 25$)  and  for  a  relatively  low  level  (T'=.03$).  As 

shown  in  the  figure,  freestream  turbulence  level  has  only  a 
slight  effect  on  transition  Reynolds  number  based  on  surface 
conditions,  (ReXt ^wsUext//vw ‘  Consistent  with  the  measure¬ 
ments,  computed  (ReXt)w  is  approximately  inversely  propor¬ 
tional  to  (Tw/Te)2;  computed  values  of  (Rext)w  generally  are 
within  about  20$  of  corresponding  measured  values. 

A  second,  more  subtle,  test  of  the  theory  is  to  apply  the 
equations  to  rough-wall  aerodynamic  ooundary  layers.  Surface 

cooling  can  actually  reverse  its  stabilizing  role  to  one  of 
destabilization  when  the  surface  is  rough.  The  explanation 
for  this  phenomenon  is  well  understood.  The  destabil  zation 
occurs  because  cooling  thins  the  boundary  layer,  thus  making 
the  roughness  look  larger  relative  to  a  boundary-layer  thick¬ 
ness  such  as  momentum  thickness,  6.  Hence,  the  surface  looks 
rougher,  whereby  transition  occurs  earlier. 

To  test  the  model's  ability  to  predict  this  phenomenon,  com¬ 
putations  were  performed  for  Mach  5  airflow  past  a  roughened 
spherical  body  with  nose  radius  of  2.5  inches;  roughness 
heights,  k,  between  1.5  mils  and  10  mils  were  used  in  the 
computations.  Freestream  unit  Reynolds  number  ranged  from 
3  to  10  million  per  foot.  (The  compressible  version  of  the 
model  equations  was  used  for  these  calculations).  The  ratio 
of  surface  temperature  to  freestream  total  temperature,  T.  , 
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was  varied  from  0.2  to  0.8.  Figure  2  compares  results  of 

the  computations  with  a  correlation  of  experimental  data'. 

Consistent  with  the  measured  effect,  momentum  thickness 

Reynolds  number  at  transition,  ReQ,  increases  with  de- 

°t 

creasing  surface  temperature . 


5  10  20 


Figure  2.  Comparison  between  computed  and  measured 
effects  of  surface  cooling  on  rough-wall 
aerodynamic  boundary-layer  transition. 
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3-2 


HEATED  HYDRODYNAMIC  BOUNDARY  LAYERS  ON  SMALL  BODIES 


The  second  round  of  applications  is  to  a  small  heated  hydro- 
dynamic  body.  The  body,  shown  in  Figure  3,  is  a  Parsons- 
Goodson-Goldschmied8  (PPG)  minimum-drag  hull  shape  which  will 
be  referred  to  as  the  H-2  body;  the  PPG  body  parameters  are 
given  in  Appendix  B.  The  objective  of  the  computations  was 
to  determine  incipient  transition  conditions  for  various 
surface  temperatures.  By  definition,  incipient;  transition 
occurs  when  the  transition  point  is  located  at  the  maximum 
body  radius. 


All  computations  were  performed  with  the  following  conditions 
specified: 


Ambient  Temperature,  T0 
Roughness  Height,  k 
Freestream  Turbulence, 
Freestream  Dissipation  Rate,  Ge 


=  71°F 
=  16  yin 
=  .01% 

=  .0185  Ug/yg 


The  pressure  distribution  was  obtained  from  the  Douglas- 
Neumann9  potential  flow  program;  Figure  4  shows  the  computed 
pressure  coefficient,  Cp,  defined  by 


■p  =  (p-p„)/(i/2pu£> 


(23) 


In  the  figure,  z  is  axial  distance  from  the  stagnation  point. 
Three  series  of  computations  were  performed  in  which  surface 
temperature  was  held  constant  at  values  of  71  f>  7*3  F>  and 
8l°F.  Incipient  transition  Reynolds  number,  Re^.  (based  on 
freestream  flow  conditions  and  V1/3  where  V  is  body  volume', 
was  computed  for  each  value  of  Tw>  Results  of  the  calculations 
are  summarized  in  Table  1;  in  addition  to  Rey,  the  table  gives 
freestream  velocity,  Uw,  transition  Reynolds  number  based  on 
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Figure  Computed  pressure  coefficient  for 

the  small  hydrodynamic  body 


In  Equation  24,  s  =  4  feet  (the  maximum  radius  point),  and 

1/ 

q  (s)  is  the  local  heat  transfer.  Figure  5  shows  computed  Rev 
w  * 

as  a  function  of  T  -T  and  q  .  As  can  be  seen  from  the  figure, 

we  w 

heating  stabilizes  the  boundary  layer  and  hence  increases  the 
value  of  Rey  This  prediction  is  qualitatively  consistent 
with  the  observed  stabilizing  effect  of  surface  heating  when 
the  fluid  is  a  liquid.  The  magnitude  of  the  effect  is  also 
realistic"1" . 


Table  1.  Computed  Incipient  Transition  Conditions  for 
the  H-2  Body  for  Varying  Surface  Temperature. 


-  ,kwatts^ 
qw('  -ft2-) 


U  (knots) 

oo 

10_6Rev 

25.3 

9.4 

35.0 

13.0 

47.4 

17.6 

10-6Re, 


18 . 2 
24.5 
33.8 


An  interesting  feature  of  the  predictions  is  revealed  in 

Figure  5-  Specifically,  Rev  increases  approximately  linearly 

with  Tw~Te  for  tne  range  of  temperatures  considered.  However, 

Re,,  increases  much  less  rapidly  with  q  .  Hence,  increasing 
V  w  2 

the  heating  rate  beyond  about  2  kwatts/ft*"  may  not  yield 

substantial  increases  in  Rey 
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3.3  HEATED  HYDRODYNAMIC  BOUNDARY  LAYERS  ON  LARGE  BODIES 

The  final  application  of  the  transition  model  is  to  a  large 
heated  hydrodynamic  body.  Again,  a  PPG  hull  shape,  referred 
to  as  .he  R-9a  body,  is  used  (see  Appendix  B).  The  objective 
of  these  computations  was  to  determine  transition  point  lo¬ 
cation  for  specified  flow  conditions. 

Two  computations  were  performed;  in  one,  the  freestream  flow 
speed,  UB,  was  35  knots  and  in  the  other  UOT  was  30  knots.  In 
both  computations,  the  following  conditions  were  imposed: 

Ambient  Temperature,  Tg 
Surface  Temperature,  Tw 
Roughness  Height,  k 
Freestream  Turbulence,  T' 

Freestream  Dissipation  Rate, 

Figure  6  shows  the  Douglas-Neumann-computed  pressure  distribu 

tion. 


*  55°F 
=  85°F 
=  16  yin 
=  .01% 

=  .0185  U2/ue 


Table  2  and  Figure  7  summarize  results  of  the  computations. 
For  both  flow  speeds,  transiton  occurs  well  upstream  of  the 
maximum  radius.  Computed  arclength  transition  Reynolds  num¬ 
bers,  Res  ,  are  not  much  larger  than  the  value  achieved  on 
the  H-2  body  with  qw  =  1.30  kwatts/ft2  (see  Table  1).  Noting 
the  heating  rates  involved,  the  predictions  are  consistent 
with  results  of  the  preceding  subsection. 
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4.  SUMMARY  AND  CONCLUSIONS 


Results  reported  in  Section  3  demonstrate  the  transition 
model’s  accuracy  for  boundary  layers  with  surface  heat 
transfer,  provided  nonuniform  viscosity  effects  are  taken 
into  account  (Subsection  2.2).  With  no  parameter  adjustment, 
the  model  accurately  predicts  transition  Reynolds  number 
for  (a)  incompressible  smooth-wall  aerodynamic  boundary 
layers  with  surface  heating  and  cooling,  (b)  cooled  rough- 
wall  boundary  layers  on  blunt  bodies  in  a  hypersonic  air- 
stream,  and  (c)  heated  hydrodynamic  boundary  layers  on  small 

bodies . 

Results  of  the  hydrodynamic  computations  agree  with  the 
measured  strong  boundary-layer  stabilization  attending 
small  amounts  of  surface  heating.  The  model  also  predicts 
that  large  amounts  of  heating  are  not  significantly  more 
effective  in  delaying  transition  than  the  smaller  rates 
considered  in  the  H-2  body  computations  of  Subsection  3.2. 
The  larger  R-9a  body  computations  support  this  prediction. 


. 


APPENDIX  A 

THERMODYNAMIC  PROPERTIES  OF  WATER 

The  following  thermodynamic  properties  of  water  are  perti¬ 
nent  to  boundary  layer  transition:  mass  density  (p;,  speci¬ 
fic  heat  (C  ),  thermal  conductivity,  (<),  molecular  viscosity 
(y),  rnd  laminar  Prandtl  number  (PrL).  Values  used  in  this 
study  for  these  quantities  have  been  obtained  from  .  chlicht- 
ing10and  are  valid  for  temperatures  ranging  from  about  40°F. 
to  110°F. 

Mass  density,  specific  heat,  and  thermal  conductivity  are 
approximately  constant  and  have  the  following  values: 


=  1.936 

lbf • sec2/ft4 

(Al) 

=  1 

Btu/lbf • °F 

(A2) 

=  0.35 

Btu/f t • hr • °F 

(A3) 

Viscosity  and  Prandtl  number,  by  contrast,  are  strongly  tem¬ 
perature  dependent  over  this  range  of  temperatures.  An  ap¬ 
proximate  polynomial  fit  to  the  Schlichting  data  was  used  to 
calculate  y  in  the  computations  of  Section  3,  i.e., 

-6  -  8 

y  =  7.943*10  (T/600)  lbfsec/ft2  (A4) 

with  temperature,  T,  in  degrees  Rankine.  Finally,  the  laminar 
Prandtl  number  follows  from  its  definition,  namely,  PrL=yCp/<, 
wherefore , 

PrT  =  2.5KT/600)"8  (A5) 

U 

Figure  A1  compares  Equations  A4  and  A5  with  the  measured  de¬ 
pendences  of  y  and  Pr^  upon  temperature. 


* 
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APPENDIX  B 

HYDRODYNAMIC  BODY  SHAPE  PARAMETERS 


The  two  hydrodynamic  bodies  analyzed  in  this  study  are 
specified  in  terms  of  the  eight  dimensionless  parameters 
used  by  Parsons,  Goodson,  and  ^cldschmied  .  The  values  of 
these  parameters  are  as  follows: 

Parameter  H-2  Body  R-9a  Body 


f 

4.450 

5.000 

r 

xm 

0.427 

0.600 

k. 

1 

2.715 

4.000 

rn 

0.260 

0.050 

ri 

0.740 

0.700 

si 

1.650 

1.750 

Xi 

0.595 

0.800 

t 

0.180 

0.150 

Additionally,  the  H-2  body  is  9-457  feet  long  and  has  a  volume 
of  14. 5  cubic  feet;  the  R-9a  body  is  55  feet  long  and  has  a 
volume  of  2500  cubic  feet. 
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